NASA Technical Memorandum 105246 



Towards the Development of Micromechanics 
Equations for Ceramic Matrix Composites 
via Fiber Substracturing 


P.L.N. Murthy and C.C. Chamis 
Lewis Research Center 
Cleveland, Ohio 


(MA ? TOWARDS I He. L^ : V f' LOPMl NT OF 1> 

MTCft jMi _r.HA.NTC S F wOAT I n N$ FLf< C t RAM IC MATRIX 
C ' J f: I T r «. V I A F 1 n SU r- STkUC T 0°. 1 **G ( *'i A S A ) 

' n C SC L 1 in unci «s 

uC/FtSi 


February 1992 




TOWARDS THE DEVELOPMENT OF MICROMECHANICS EQUATIONS FOR 
CERAMIC MATRIX COMPOSITES VIA FIBER SUBSTRUCTURING 

P.L.N. Murthy and C.C. Chamis 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

ABSTRACT 

A generic unit cell model which includes a unique fiber substructuring concept is proposed 
for the development of micromechanics equations for continuous fiber reinforced ceramic 
composites. The unit cell consists of three constituents: fiber, matrix and an interphase. In the 
present approach, the unit cell is further subdivided into several slices and the equations of 
micromechanics are derived for each slice. These are subsequently integrated to obtain ply level 
properties. A stand-alone computer code containing the micromechanics model as a module is 
currently being developed specifically for the analysis of ceramic matrix composites. Towards 
this development, equivalent ply property results for a SiC (silicon carbide fiber)/Ti-15-3 
(Titanium matrix) composite with a 0.5 fiber volume ratio are presented and compared with 
those obtained from customary micromechanics models to illustrate the concept. Also, 
comparisons with limited experimental data for the ceramic matrix composite, SiC/RBSN 
(Reaction Bonded Silicon Nitride) with a 0.3 fiber volume ratio are given to validate the 
concepts. 

SYMBOLS 

C heat capacity 

d f fiber diameter 

E,G Young’s (Normal) and shear moduli 

f,m,i subscripts for fiber, matrix, and interphase 

K thermal conductivity 


k 


volume fraction 


t subscript for ply or slice (subply) 

n Poisson’s ratio 

T,AT melting temperature, temperature difference 

1,2,3 subscripts for material axis. 1 along the fibers 2 and 3 are transverse to the fiber 

a coefficient of thermal expansion 

e strain 

p density 

a stress 

INTRODUCTION 

Composite micromechanics is an extensively explored field where the primary thrust is to 
evaluate the behavior of the composite as a function of the behavior and interaction of the 
constituents (fiber and matrix). There are two basic approaches to the micromechanics of 
composite materials: (1) mechanics of materials and (2) elasticity. The principles and the 
assumptions behind these approaches are well known and are standard topics of many text books 
[Ref. 1, for example]. 

In-house research over the past twenty years has focused on the mechanics of materials 
approach and has culminated into several computer codes for composite micromechanics and 
macromechanics. Notable among these are ICAN (Integrated Composite ANalyzer) and 
METCAN (METtal Matrix Composite ANalyzer). ICAN [Ref. 2] was primarily developed for 
predicting the behavior of polymer-matrix composites. The micromechanics of ICAN are based 
upon a unit cell consisting of the fiber and matrix that is arranged in a regular square 
array pattern (Fig. 1). A recently developed code using the same philosophy is METCAN 
[Ref. 3] which was tailored for analyzing micromechanical as well as macromechanical behavior 
of metal matrix composites. The unit cell model for METCAN is slightly different in that it 
consists of three distinct constituents fiber, matrix and an interphase. The third constituent, the 
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interphase is unique to this class of materials which can develop due to a possible chemical 
reaction between the fiber and matrix. 

Currently, a unique fiber substructuring concept is being explored to further describe the 
micromechanics of fiber composites with finer local detail. In this approach the customary unit 
cell model is substructured into several slices (Fig. 2). Each slice is treated as the smallest 
representative unit and micromechanics equations are derived for each slice by applying the well 
known principles of the mechanics of materials approach. The basic philosophy is applicable to 
any type of continuous-fiber-reinforced matrix composite. However, the current emphasis is on 
the development of a dedicated stand-alone computer code using the above mentioned concepts 
for predicting the behavior of ceramic matrix composites (CMCs). The code is thus named 
CEMCAN for CEramic Matrix Composite ANalyzer. 

The current approach has several advantages over other conventional methods. A few of 
these advantages provide the analyst/designer the capability to better distribute local stresses, 
the ability to specify various degrees of bond around the fiber circumference, and to account for 
the fiber breaks and matrix cracking. These are some of the issues that must be addressed to 
adequately describe the micromechanical behavior of ceramic matrix composites. 

The objective of the present paper is twofold. First, to describe the fiber substructuring 
concepts and the micromechanics equations that are embedded in CEMCAN; and second, to 
describe the computer code CEMCAN, its current features and capabilities and some illustrative 
examples to demonstrate the versatility of the code. 

FIBER SUBSTRUCTURING AND MICROMECHANICS 

The primary objective of composite micromechanics is to determine the equivalent elastic 
moduli of a composite material in terms of the elastic moduli of the constituent materials. An 
additional and complementary objective of the micromechanics approaches is to determine the 
strengths of the composite material in terms of the strengths of the constituent materials. Other 
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properties of interest are the composite thermal expansion coefficients, thermal conductivities 
and heat capacity. 

The derivation of equivalent composite ply properties using micromechanics theory starts 
with the identification of a representative volume element or unit cell. It is the smallest region 
or piece of material over which the stresses and strains are macroscopically uniform. However, 
within the unit cell the stresses and strains are nonuniform due to the heterogeneity of the 
material. The unit cell can consist of a fiber, matrix and/or an interphase. It is assumed that 
these cells are arranged in a regular pattern. Furthermore, it is assumed that the unit cell be 
either a square or a hexagon depending upon the chosen array pattern. Equivalent properties for 
the ply are then derived in terms of the constituent material properties based on the mechanics 
of materials approach. Other assumptions involved in this approach are: (1) fiber and matrix 
are subjected to the same strain in the fiber direction of a unidirectional fibrous composite, and 
(2) the same transverse stress is applied to both the fiber and matrix in the direction transverse 
to the fiber. 

The details of the representative volume element chosen previously in the development of 
ICAN [Ref. 2] and METCAN [Ref. 3-5] are shown in Fig. 1. The unit cell for ICAN consists of 
two distinct regions A (matrix only) and B (fiber and matrix). The unit cell for METCAN, 
however, may consist of a maximum of three regions, A, B, and C. The region A consists of 
matrix only, the region B consists of a combination of matrix and interphase, and the region C 
consists of a combination of all three constituents: fiber, matrix and interphase. Note, that in 
the following discussion the interphase is treated as a separate constituent with distinct 
properties. Thus, it can either represent a zone formed due to a chemical reaction between the 
fiber and matrix or a separate layer provided intentionally to prevent such a reaction. The 
different regions facilitate the representation of the nonuniformity in the local stress distribution. 
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In the application of this approach to CMCs the same type of array pattern is assumed. 
However, the unit cell is further subdivided into several slices and the equations of 
micromechanics are derived for each slice. These are subsequently integrated to obtain unit 
cell/ply level properties. The modeling details of the substructured unit cell are illustrated in 
Fig. 2. To derive the equivalent slice properties consider a slice taken out from the unit cell. 
Let d { , d m , and d } be the fiber, matrix, and interphase widths. Let h be the height of the 
slice and s be the total width of the slice. Since it is customary to express the equivalent 
properties in terms of the properties of each of the constituents and their respective volume 
ratios, let k f , k m , and k. be the fiber, matrix and the interphase volume ratios respectively. 
Then by definition 


d <. k _ 2 V k . _ 2d i 

— * A m 1 


(i) 


s s s 

Now that the preliminary parameters have been established; the mechanics of materials 
approach is applied to generate the composite mechanical and thermal properties. 

Mechanical Properties 

By applying the force equilibrium in the longitudinal direction the following equation can 
be written. 


<r |U hs = (2<r mll d m + 2ff in dj + <r fll d f )h (2) 

Also, based on the assumptions for the mechanics of the materials approach one can write the 
following equation. 


g lll - e mll = e ill - e fll ( 3 ) 

With the aid of Eqs. (1) to (3) the longitudinal modulus, E^, can be written as 

Ejn = k f E fll + k m ®mll + ^i^iii ( 4 ) 

The Poisson’s ratio in 1-2 direction can be derived by considering the strain in the transverse 
direction (22) as the sum of the individual strains in the constituents. Accordingly, 


5 


c 122 s = 2e m22 d m + 2 ®i22 d i + c f22 d f 


( 5 ) 


Also, by definition 


i/; 


jl2 


®j 22 
E jll 


( 6 ) 


where j can represent any of the subscripts t, m, i, or f. 

The combination of Eqs. (1), (5), and (6) leads to the following equation for the major 
Poisson’s ratio 


I/ I12 = + ^m I/ ml2 + ^i I/ il2 


( 7 ) 


The modulus in the transverse (2-2) direction can be derived by assuming that the same 
transverse stress, a is applied to the fiber, matrix and interphase. 

ff (22 = ff m22 = a i22 ~ a i22 
With the aid of Eqs. (5) and (8) and noting that 


'j 22 


ff j22 

Ej22 


( 9 ) 


one can derive the following equation for the transverse modulus in the direction 2-2. 


E 


E m 22Ei22 E f22 


122 


km E i22 E f22 + kj E f22E m 22 + kfE m 22Ei22 

The in-plane shear modulus in (1-2) direction is determined by assuming that the shearing 

stresses on the fiber and the matrix are the same. In this respect the derivation of G {12 is 

similar to that for E^. The equation for G^ 2 can be shown as 

G ml2 G il2 G fl2 


( 10 ) 


G 


(12 


km G il2 G fl2 + ki G fl2 G ml2 + ^f G ml2 G il2 

The remaining mechanical properties can be derived in a similar fashion and are given by 


( 11 ) 


E m 33 E i 33 Ef 33 


-(33 


kmEi33Ef 33 + k i E f33 E m 33 + k f E m33 E i33 


( 12 ) 
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* 

G ml3 G il3 G fl3 

(13) 

^*113 

k m G il3 G fl3 + 

k i G fl3 G ml3 

+ k f G ml3 G il3 

n - 

< 

3 m23 G i23 G f23 

(14) 

^*123 ~ 

k m G i23 G fZ3 + 

k i G f23 G m23 

+ k f G m23 G i23 


"f 13 = k f"fl3 

+ k m"ml3 + 

k i"i 13 

(15) 


"|32 = k f"f 32 

+ k m"m32 + 

k i"i 32 

(16) 


The Poisson’s ratios in 2-3, 2-1, and 3-1 directions are given by the reciprocity relations 


"123 = 


"|32 E |22 

E «33 


(17) 


"121 - 


"| 12 E |22 


E 


111 


(18) 


"*31 “ 


"|13 E |33 

E lll 


(19) 


Thermal Properties 

The equivalent thermal conductivity in the (l-l) direction can be derived by equating the 
total heat flow rate across the cross section to the sum of the individual rates in each 
constituent. The heat flow rate can be calculated using the Fourier’s law of heat conduction. 
By following the above steps the heat flow rate equilibrium equation can be written as 


K m sh£I - k tu d,K^I * <■ K 

dL dL dL 


i > dT 

‘“ d ' h dL 


( 20 ) 


where (dT/dL) is the thermal gradient across a typical length of the cell L in the (1-1) 
direction (Fig. 3). Simplification of the above equation leads to the following; 
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( 21 ) 


K,n = k f K fn + k m K mll + kjKj X1 
The equivalent thermal conductivities in the (2-2) can be derived by noting that the total heat 
flow rate across each constituent in the (2-2) direction should be same. Thus from Fig. 3 the 
following can be written: 


Q = IWLh) 


- Km(Lh) 


T 5 - To 


t 3 -t 2 


- K m22 (Lh) 


= K i22 (Lh) 


t 5 -t 4 


m 


T 2 -Tj 


= K i 22 (Lh) 


t 4 -t 3 


d; 


= K m 22 (Lh) 


d; 


Tx - T 0 


( 22 ) 


in 


where the T Q through T 5 are l.he temperatures at various cross sections along the slice as 
shown in the Fig. 3. The above equation can be rearranged as 


K| 22 “ K m22 


T5 - T 4 


T 5 - To 


m 


T 5 -t 4 


(23) 


m 


[(t 5 - T 4 ) + (t 4 - T 3 ) + (t 3 - T 2 ) + (t 2 - T>) + (T x - T 0 )J 

By observing the fact that (T 4 - T 3 ) = (T 2 - Tj) and (T 5 - T 4 ) = (T x - T 0 ), the above can be 
further simplified into 


K 


(22 


2+2 

[T 4 - T 3 j 

+ 

t 3 - t 2 | 


[t 5 - t 4 J 


i T S - T 4jJ 


K 


m22 


(24) 


With the aid of Eqs. (1) and (22), Eq. (24) can be written as K {22 = 
which leads to 


K 


m22 


kjKm22 k f K m22 


k m + 


K; 


22 


K 


f 22 


K 


(22 


^m22^i22^f22 


k m K i22^f22 + k i^f22^m22 + k f^m22^i22 


(25) 
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The thermal conductivity along the (3-3) direction can be derived in a similar manner and is 
given by 


K 


^m33^i33^f33 


(33 


^m^i33^f33 + ^i^f33^m33 + ^f^m33^i33 
The equivalent heat capacity Cj can be derived by equating the sum of the heat contained in 

the individual regions to the heat contained in the slice based on the equivalent heat capacity: 


( 26 ) 






kjPjCi + k f p f C f 


(27) 


Pi 

where p ^ is the equivalent density of the slice and p m , p., and p f are respectively the 
densities of the matrix, interphase and fiber. The equivalent density can be expressed in terms 
of the individual densities by the simple rule of mixtures: 


P( = Vf + k m Pm + k i^i ^ 28 ) 

The longitudinal thermal expansion coefficient can be derived by noting that the sum of the 
forces in the longitudinal direction should be zero: 


( 2d m®mll + 2d i®ill + d f ll)°lll A T 

= "( 2d m a mll®mll + 2d i a ill®ill + d f a f ll^fll) = 0 


(29) 


where AT is the temperature difference. Simplification of Eq. (29) leads to the following 
expression for the longitudinal thermal expansion coefficient: 

k m a mll^mll + k i a ill^ill + k f a fll^fll 


"ill = 


E 


(11 


(30) 


The transverse thermal expansion coefficient can be derived by equating the total strain in the 
(2-2) direction to the sum of the individual strains in each region: 

<*|22 s AT " 2d m"m22 AT + 2d i"i22 AT + d f"f22 AT ( 31 ) 

With the aid of Eqs. (l) and (31) a^ 2 can be written as 

"(22 = k f"f22 + k m"m22 + k i"i22 ( 32 ) 
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The expression for a^ 3 is similar to the above and is given by 

°«33 = k f a f33 + k m a m33 + k i a i33 ( 33 ) 

This set of equations define the equivalent properties for a typical slice with three regions. The 
procedure to obtain the equivalent properties of the representative volume element (unit cell) is 
analogous to that of obtaining the properties of a typical laminate from the ply properties. Here 
each slice is treated as a subply within the unit cell. 

CEMCAN COMPUTER CODE 

These new micromechanics equations have been programmed as a separate 
micromechanics module and integrated into a computer code similar to that of the earlier ICAN 
and METCAN codes. The new code is named CEMCAN and its flow chart is shown in Fig. 4. 
Not all aspects of the code as shown in the figure, are currently operational; a brief description 
follows. 

The material behavior model depicted in the center describes the constituent property as a 
function of several independent variables such as the current temperature, stress and time. On 
the basis of the current properties, the equivalent properties for the slice are calculated using the 
micromechanics equations. The left side of the Fig. 4 shows the progressive integration or 
synthesis of the various stages in CEMCAN starting from the slice level properties and ending 
with the description of the laminate level properties. The details involved are described briefly 
in the following paragraphs. Some of these details are of text book nature and therefore the 
reader is advised to consult text books on composite mechanics [Ref. 1, for example] for a 
complete in depth discussion. 

Step 1. Properties of a lamina from the slice level properties . 

This involves integration of slice level properties through the unit cell to obtain the 
equivalent properties of the unit cell. Since it is the smallest representative volume element, 
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these properties can be considered as equivalent properties of a lamina. The integration from a 
slice to unit cell is achieved by using the classical laminate theory. 

Step 2. Properties of a single ply from the lamina properties . 

To obtain ply properties from lamina properties two cases must be considered. First, if 
there is a single fiber through the thickness of the ply, then the ply properties are identical to 

i 

those of the lamina. However, if there are a number of fibers through the thickness of a single 

I 

ply, then this ply is treated as a laminate made of a number of single fiber laminae. 

Consequently, the ply properties can be obtained by using once again the classical laminate 

theory. i 

Step 3. Properties of a laminate made of several plies . 

The next step of integration involves the determination of global (laminate or structural 
level) properties from the single ply properties using macromechanics. 

Step 4. Global or Laminate response to loads . 

The global or laminate response at a point through the thickness consists of the strains 
and curvatures developed due to applied mechanical and thermal loads. This is accomplished 
using the well known principles of classical laminate theory. 

The right hand side of the figure (Fig. 3) shows the successive reduction or decomposition 
of the laminate level response (strains and curvatures) to the slice level response. This involves 
going through the same steps as above in a reverse order. 

Currently the micromechanics equations for the slice and integration to the ply level 
properties as well as decomposition of the ply response to the slice level response are completely 
programmed and are operational. In addition, the code can accommodate circular and 
arbitrarily shaped fibers by altering slicing in the horizontal and vertical directions (Fig. 5). 

Additional features including the consideration of material nonlinearities and constituent 
response within the slice are planned for the future. 
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TYPICAL RESULTS AND DISCUSSION 

The first set of results ensue from a verification exercise in which the proper formulation 
of the equations is being checked using the established micromechanics in ICAN and METCAN. 
To verify /illustrate the code’s features and capabilities a typical composite system SiC (Silicon 
carbide fiber)/Ti-15-3 (Titanium alloy matrix) is chosen. The reasons for this particular choice 
are (1) the selected composite system is of current interest and (2) at this point the intention is 
only to verify the micromechanics equations developed using the new fiber substructuring 
concept. The typical properties for these constituents based upon the data given in Ref. [6] are 
shown in Table I. Keeping a constant fiber volume ratio of 0.5, composite properties for the 
SiC/Ti-15 are generated using CEMCAN as well as ICAN and METCAN. The resulting 
predictions for composite moduli, Poisson’s ratios, thermal conductivities and the coefficients of 
thermal expansion are shown in Figs. 5 to 8. 

Mechanical Properties 

Unidirectional ply normal and shear moduli for the SiC/Til5 composite system (with a 
fiber volume ratio (FVR) of 0.5) are shown in Fig. 6. All three codes predicted almost the same 
value for the ply longitudinal modulus. There appears to be much better agreement between 
CEMCAN and METCAN for the remaining moduli. ICAN consistently over predicted with the 
exception of the through the thickness shear modulus G j 23 . This is to be expected because, the 
micromechanics equations embedded in ICAN are strictly for polymeric matrix composites which 
involved certain simplifying assumptions [Ref. 7]. Similar trends are noticed in Fig. 7 also, where 
the predictions of the Poisson’s ratios for the same composite system are compared. However, 
the differences, if any, are minimal. 

Thermal Properties 

Unidirectional ply thermal expansion coefficients for the composite system examined appear in 
Fig. 8. With the exception of the longitudinal coefficient of thermal expansion a, all three 
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codes are in good agreement. The longitudinal thermal expansion coefficient predicted by 
METCAN is significantly higher (18 percent). The thermal conductivities predicted by all three 
codes are in excellent agreement as shown in Fig. 9. 

Mechanical Response 

The stress variation in a ply subjected to combined in plane loading is shown in Fig. 10. 
The benefits of fiber substructuring becomes clear from this example. The conventional 
micromechanics model in ICAN predicted a constant stress in the ply. In contrast, CEMCAN 
predicts a smooth variation of stresses through the ply thickness. 

Thermal Response 

The variation in the in plane stresses due to a 100 °F (311 K) increase in temperature is 
shown in Fig. 11. Once again, due to the fiber substructuring, CEMCAN is able to predict 
smooth variations in the in-plane stresses as opposed to a conventional model which predicts no 
stresses at all. 

Verification of the CMC Micromechanics 

The experimental data for a 0.3 FVR SiC/RBSN ceramic matrix composite system 
reported in Ref. [7] is used to conduct a limited validation study. More validation studies will 
be provided in future reports on the subject. Table II lists the properties of the constituents 
used for this study. The CEMCAN predictions are compared to the experimental data as shown 
in Fig. 12. Note that the properties are normalized with respect to the theoretically predicted 
value. As seen in the figure the longitudinal modulus as well as the major Poisson’s ratio are in 
excellent agreement with the experimentally determined values. This is not very surprising 
because of the fiber dominated behavior in the longitudinal direction. The transverse modulus, 
the shear modulus and the minor Poisson’s ratio are all over predicted. The experimentally 
observed values are about 50 percent or less than the theoretical predictions. Once again this 
should be expected because the theoretical model considered perfect bond between the fiber and 
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matrix. In reality for CMC materials a very weak or nonexistent bond is quite common. To 
test this premise, a case was simulated utilizing an interphase with a thickness 5 percent of the 
diameter of the fiber and with negligible properties (listed in Table II). The results are shown in 
Fig. 11. The predicted properties now represent a lower bound for the composite. As seen from 
the figure, the experimentally measured values for the transverse direction clearly fall between 
the perfect bond and weak bond simulations. The experimental values are however much closer 
to the simulated results with weak bond. With the present model one can also try partial 
bonding as well as a host of other combinations. 

SUMMARY AND FUTURE PLANS 

Based on the results presented in this study, it can be concluded that: 

1. The prediction of ply mechanical and thermal properties agreed very well with the 
existing models in ICAN and METCAN; lending credence to the fiber substructuring approach. 

2. The verification case presented for SiC/RBSN composite system showed good 
agreement with the experimental results. Also, methods to establish upper and lower bounds to 
the mechanical properties are illustrated. 

3. Fiber substructuring can capture greater local detail than conventional unit cell based 
micromechanical theories. Therefore, it offers promise in simulating complex aspects of 
micromechanics in ceramic matrix composites (CMC) such as (a) various degrees of bond around 
the fiber circumference and along the length, (b) fiber breaks and matrix cracking, (c) effects of 
these on ply thermal and mechanical properties/response. 

4. The approach is general and versatile. It can be applied to any type of continuous fiber 
reinforced matrix composite. 

The response decomposition from the slice (subply) level to the local (in each constituent) 
level is currently under development. This will enable predictions of local stresses (fiber, matrix, 
and interphase regions) within the slice. A nonlinear material behavior model to describe the 
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current constituent properties as a function of the reference properties and the current response 
is also under development. 
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Table I. Properties Chosen in the Study for SiC Fiber and Ti-15 Matrix 


SiC fiber properties 


Ti-15 matrix properties 


E f 

62.0 Mpsi 

427 GPa 

E m 

12.3 Mpsi 

85 GPa 

G f 

23.8 Mpsi 

164 GPa 

G m 

4.7 Mpsi 

32 GPa 


0.3 



0.32 


Kf 

0.75 Btu/hr/in./°F 

1.3 w.mK 

K f 

0.39 Bt,u/hr/in./°F 

0.67 w/mK 


1.8 ppm/°F 

3.2 ppm/K 

a_ 

m 

4.5 ppm/°F 

8.1 ppm/K 

c f 

0.29 Btu/lb 

1.2 kJ/kg K 

G m 

0.12 Btu/lb 

0.5 kJ/kg K 

P[ 

0.11 lb/in. 3 

3 Mg/m 3 


0.17 lb/in. 3 

4.8 Mg/m 3 

T f 

4870 °F 

2961 K 

T 

m 

1800 °F 

1255 K 

<L 

0.0056 in. (0.14 mm) 






Table II. Properties Chosen in the Study for SiC/RBSN Composite System 


SiC fiber properties 


RBSN matrix properties 


Interphase 


E f 

56.6 Mpsi 

390 GPa 

E m 

15.9 Mpsi 

110 GPa 

E i 

ioo psi 

0.7 GPa 

G f 

24.2 Mpsi 
0.17 

0.0056 in. (0.14 mm) 

150 GPa 

G m 

6.5 Mpsi 
0.22 

45 GPa 

G i 

•'i 

40 psi 
0.22 

0.27 GPa 


M 

111 

1 

m 

i: 

M 

in 

B 

111 


Square Array 



ICAN Unit Ceil METCAN Unit Cell 

Figure 1. — Unit cell square array concepts In 
ICAN and METCAN codes. 
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Unit cell with fiber, matrix, and 
interphase divided into N|_ slices 
(subplies) 


hV 
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L 

rL 

j 

T 





1 

u 
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► 


Individual slice showing three regions: 
fiber, matrix, and interphase 



(1-1) direction 



(2-2) direction 


Figure 2. — Ply/fiber substructuring concept for Figure 3. — Modeling details for the slice for deriving the equiv- 

ceramtc matrix composite micromechanics - alent thermal conductivities, 

mechanics of matenals approach. 


Laminate Laminate 



Slice 


Figure 4.— Integrated analysis approach embedded in CEMCAN. 
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(Note: 1 1 Is along the fiber direction) 


Figure 5. — Different fiber shapes and slicing arrangements 
available in CEMCAN. 
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(0.5 fiber volume ratio) 
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bb ican 

CEMCAN 
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Figure 6.— Comparison of ICAN, CEMCAN, and M ETC AN 
predictions of Young's moduli and shear moduli for unidirec- 
tional SiC/Til 5-3 composite at .5 fvr. 


Figure 7.— Comparison of ICAN, CEMCAN, and M ETC AN 
predictions of Poisson's ratios for a unidirectional SiC/Ti15-3 
composite at .5 fvr. 
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Figure 8. — Comparison of ICAN, CEMCAN, and METCAN 
predictions of thermal expansion coefficients for a unidirec- 
tional SiC/Til 5-3 at .5 fvr. 


Figure 9. — Comparison of ICAN, CEMCAN, and METCAN 
predictions of thermal conductivities for a unidirectional 
SiC/Til 5-3 composite at .5 fvr. 
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Figure 10.— Variation of longitudinal stress within the ply as 
predicted by CEMCAN for a unidirectional SiC/Ti15-3 
composite at .5 fvr under a combined in~plane loading. 


Figure 11.— Variation of the in-plane stresses within the ply as 
predicted by CEMCAN for a unidirectional SiC/Ti15-3 com- 
posite at .5 fvr under a 1 00 °F uniform heating. 
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Figure 1 2.— Mechanical properties of SiC/RBSN composite 
(0.3 fvr): CEMCAN vs experiments. 
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